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A simple and efficient approximation scheme to study electronic transport characteristics of 
strongly correlated nano devices, molecular junctions or heterostructures out of equilibrium is pro¬ 
vided by steady-state cluster perturbation theory. In this work, we improve the starting point of 
this perturbative, nonequilibrium Green’s function based method. Specifically, we employ an im¬ 
proved unperturbed (so-called reference) state p s , constructed as the steady-state of a quantum 
master equation within the Born-Markov approximation. This resulting hybrid method inherits 
beneficial aspects of both, the quantum master equation as well as the nonequilibrium Green’s 
function technique. We benchmark the new scheme on two experimentally relevant systems in the 
single-electron transistor regime: An electron-electron interaction based quantum diode and a triple 
quantum dot ring junction, which both feature negative differential conductance. The results of the 
new method improve significantly with respect to the plain quantum master equation treatment at 
modest additional computational cost. 

PACS numbers: 71.15.-m, 71.27+a, 73.63.-b, 73.63.Kv 


I. INTRODUCTION 

Electronic transport in the realm of molecular scale 
junctions and devices has become a subject of intense 
study in recent years. 1-7 Nowadays the controlled as¬ 
sembly of structures^ via electro migration^— the con¬ 
tacting in mechanical break-junction setups^ - — elec¬ 
tronic gatin g 17 : 19,22 and measurement via scanning tun¬ 
nelling microscopy^r— have become established tools, 
ultimately opening routes from elementary understand¬ 
ing to device engineering. Prompted by these formidable 
advances in experimental techniques, the characteriza¬ 
tion of transport through e.g. molecules bound by 
anchor groups to metal electrodes/ 11,2 ' -28 heterostruc¬ 
ture s 29 ! 30 or nano structures on two- dimensional sub¬ 
strate s 29 ! 31 — has become feasible. These constitute the 
foundation for future applications in electronic devices 
based on single electron tunnelling^ quantum interfer¬ 
ence effectsspin contro l 45 ! 46 or even quantum many- 
body effects^ ! 10 ! 12 ! 47 ! 48 like Kondo^ behaviour— T— 

Typically the electronic transport through such de¬ 
vices is significantly influenced by electronic correlation 
effects, which may become large due to the reduced ef¬ 
fective dimensionality and/or confined geometries. This 
is reflected, for instance, in major discrepancies between 
experimental and theoretical current-voltage characteris¬ 
tics obtained with (uncorrelated) nonequilibrium Green’s 
function^ - — calculations based on ab-initio density 
functional theory statesThe inclusion of many- 
body effects in the theoretical description of fermionic 
systems out of equilibriu m 55 1 63 ~ — is challenging and an 
active area of current research—T— Suitable approxima¬ 
tions need to be devised in order to solve a finite strongly 
correlated quantum many-body problem out of equilib¬ 
rium coupled to an infinite environment. Typically, the 
nonequilibrium setup consists of a correlated central re¬ 
gion (system) attached to two leads (environment). 


A well-established method for treating such open quan¬ 
tum systems is by means of quantum master equations 
(Qme).— — Herein, the environment-degrees of freedom 
are integrated out and usually incorporated in a pertur¬ 
bative manner. The Qme approach allows a detailed in¬ 
vestigation of transport phenomen a 45 ’ 46 and recent self- 
consistent extensions attempt to cure some of its long¬ 
standing limitations— 

In the framework of nonequilibrium Green’s functions 
(NEGF) various schemes exist to approximately calcu¬ 
late the electronic self-energy of the correlated region, 
see e.g. Ref. IfisHs*! [87l In cluster approaches, such as 
cluster perturbation theory (CPT) and its improvement, 
the variational cluster approach (VCA)^ the whole sys¬ 
tem is partitioned into parts which can be treated ex¬ 
actly and determine the self-energy. Originally devised 
for strongly correlated systems in equilibrium, 89- — both 
approaches have recently been extended to nonequilib¬ 
rium situations in the time dependent cas o 91 ' 92 as well 
as in the steady-state. 92,93 In previous work we applied 
the steady-state CPT (stsCPT) to obtain transport char¬ 
acteristics of heterostructures ^ quantum dots^^r— and 
molecular junction s 97 ! 98 and obtained good results even 
in the challenging Kondo regime . 49 ' 94 ! 95 

A key issue in the CPT approach is to identify an ap¬ 
propriate many-body state for the disconnected corre¬ 
lated cluster in the central region, as a starting point of 
perturbation theory, the so-called reference state. Up to 
now, a common choice in stsCPT is to use an equilibrium 
state at some temperature Tg (often T$ = 0) and chemi¬ 
cal potential /rs in-between the values of the leads. Such 
an ad-hoc choice is clearly unsatisfactory. Furthermore, 
it fails to describe certain quantum interference effects 
in transport phenomena as for example so-called current 
blocking effects : 45 ! 46 ' 97 

The purpose of the present work is to improve on 
stsCPT by constructing a consistent and conceptually 
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more appropriate reference state, given by the steady- 
state reduced many-body density matrix p s obtained 
from a Qme in the Born-Markov approximation. Within 
this quantum master equation based stsCPT (meCPT), 
the ambiguity in defining ps and Ts for the central re¬ 
gion is resolved. The equilibrium case, in which ps and 
Tg coincide with those of the environment, is automati¬ 
cally included. In contrast to standard Qme approaches, 
lead induced level-broadening effects are accounted for 
and the noninteracting limit is reproduced exactly, as in 
the original stsCPT. In addition, meCPT is able to cap¬ 
ture the previously mentioned current blocking effects, 
as shown below. 

Other NEGF/Qme hybrid methods exist in the 
literaturei 67 ’ 99 ! 100 For instance, in a recent wor k 101 ’ 102 
we have proposed a so-called auxiliary master equation 
(AME) approach, whereby a Lindblad equation is intro¬ 
duced which models the leads by a small number of bath 
sites plus Markovian environments. The AME is suited 
to address steady-state properties of single impurity 
problems as encountered in the framework of nonequilib¬ 
rium dynamical mean field theory . 72 87 ] 101109 ■A® In con¬ 
trast, the meCPT presented in this work is more appro¬ 
priate to treat non-local self-energy effects which cannot 
be captured by single-site DMFT. 

This paper is organized as follows: After defining the 
model Hamiltonian in Sec. eh the meCPT is introduced 
in detail in Sec. m We present results obtained with 
the improved method for two experimentally realizable 
devices: i) In Sec. IIV A1 an electron-electron interaction 
based quantum diode, ii) and in Sec. IIV B1 a triple quan¬ 
tum dot ring junction which both feature negative differ¬ 
ential conductance (NDC). 

For ring systems, extensive Qme results and an ex¬ 
planation of the NDC in terms of quantum interference 
mediated blocking are available in Ref. I45ll46l . 


II. MODEL 

We consider a model of spin-A fermions, having in mind 
the electronic degrees of freedom of a contacted nano 
structure, heterostructure or a molecular junction. The 
Hamiltonian consists of three parts: 

n=U s + U E +'H SE . (l a ) 

i) The “system” T~L S represents the interacting central 
region i.e. the nano device or molecule consisting of 
single-particle as well as interaction many-body terms. 
It is described by electronic annihilation/creation opera- 
tors fia/fl 7 at site * = I 1 ’ ■ • ■ > N s] where N s is typically 
small and spin a = {j*, 4,},A21i We will specify the par¬ 
ticular form of i-L s in the respective results section, ii) 
The “environment” Hamiltonian % E describes the two 


noninteracting electronic leads 
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= Xm e \ka C \kc C \kv . ( lb ) 

A—1 k<7 

where c xk( 7 /c\ kcr denote the fermion operators of the in¬ 
finite size lead A with energies e\ kc , and electronic den¬ 
sity of states (DOS) p\ a {w) = jA J2 $ (w - £\k*) where 

k 

N\ —> oo are the number of levels in the leads. The 
disconnected leads are held at constant temperatures T\ 
and chemical potentials p\ so that the particles obey 
the Fermi-Dirac distribution T\, / ta )' 106 ' 107 hi) Fi¬ 

nally the system and the environment are coupled by the 
single-particle hopping 

^se = (t'xikaflc Xka ) + h.c .. (lc) 

A—1 ika 


III. MASTER EQUATION BASED CLUSTER 
PERTURBATION THEORY 

Our goal is to obtain the steady-state transport charac¬ 
teristics of the Hamiltonian "H, Eq. (Hal) in a nonequilib¬ 
rium situation induced by environment parameters like 
a bias voltage Vb or temperature gradient AT. The 
important step consists in evaluating the steady-state 
single-particle Green’s function in Keldysh space G in 
the well established Keldysh-Schwinger nonequilibrium 
Green’s function formalism ii2§riA£ In general H is both 
interacting and of infinite spatial extent. Therefore ex¬ 
plicit evaluation of G is prohibitive in all but the most 
simple cases which motivates the introduction of approx¬ 
imate schemes. 

One such scheme is CPT ) 89 i 90 in which one performs 
an expansion in a ’small’ single-particle perturbation, 
for example the system-environment coupling T~L SE of 
Eq. (flcl) . The unperturbed Hamiltonian ii s + H E can 
be solved exactly. While in the noninteracting case CPT 
becomes exact, results obtained in the presence of in¬ 
teraction are approximate and depend on the reference 
state for the unperturbed system. A common practice 
within stsCPT— is to use a pure state given by the 
equilibrium ground state |'Fo)s of the disconnected in¬ 
teracting system Hamiltonian ~H S . In a nonequilibrium 
situation, this is still ambiguous, as it depends on an ar¬ 
bitrary choice of the chemical potential ps and/or tem¬ 
perature Ts for the interacting finite system. 

The goal of this work is to provide an unambiguous 
and conceptually more rigorous criterion for the choice 
of the reference state for the interacting central region. 
Ideally, the reference state is selected such that it resem¬ 
bles best the situation of the coupled system, i.e. for the 
full Hamiltonian, Eq. (flal) in the steady-state. An appro¬ 
priate choice in equilibrium is to use the grand-canonical 
density operatoriSS p E c as reference state. In this case, 












3 


Ts and ps are uniquely determined by the equilibrium 
situation. Equivalently, pf c is given by the steady-state 
solution of a Qme in the Born-Markov approximation 
(see Sec. lIIIBl) . when coupling the system to one thermal 
environment. From this viewpoint a natural extension to 
the nonequilibrium situation is to make use of a Qme as 
well in order to obtain a consistent reference state, given 
then by the steady-state reduced density operator of the 
system p s . In this work, a second order Born-Markov 
Qme is employed, which yields the correct zeroth order 
reduced density operator p s (adjusted to 'H SE )} 11A12 
Subsequently, TL SE is included within the CPT approx¬ 
imation, 89,90 in order to obtain improved results for the 
Green’s function and in turn for the transport observ¬ 
ables. 

In summary, the meCPT method consists of the fol¬ 
lowing three main steps, analogous to a standard CPT 
treatment: 


1. Decompose the whole system into a small inter¬ 
acting central region (system) and noninteracting 
leads of infinite size (environment), see j-L s and 1 ~L E 
in Eq. (flal) . 

2. The new step introduced in this work is to solve a 
Qme for the system in order to obtain the reduced 
density operator p s , which serves as a reference 
state to calculate the cluster (retarded) Green’s 
function^ 



—*0(r)tr 



fMJja 



( 2 ) 


3. Reintroduce the system-environment coupling 'H SE 
perturbatively, see Sec. IIII Al and Eq. fl]), to deter¬ 
mine the Green’s function of the coupled system. 


A. Steady-state cluster perturbation theory 

Here we briefly recall the main, well-established CPT 
concepts and equations, as this is the starting point for 
the formalism presented in this work. For an in depth 
discussion of CPTii^ and its nonequilibrium extension 
we refer the reader to the literature ; 91 ' 93 ! 95 ! 97 

The central element of stsCPT is the steady-state 
single-particle Green’s function in Keldysh spaced 

5-CT£). w 

where R denotes the retarded, A the advanced, and K the 
Keldysh component. In the present formalism, G R ^ A ^ K 
become matrices in the space of cluster sites and depend 
on one energy variable u since time translational invari¬ 
ance applies in the steady-state. 

As explained above, in order to compute G(co) within 
stsCPT one partitions H, Eq. (flal) in real space, into in¬ 
dividually exactly solvable parts, in this case, the system 


T~t s and the environment T~L E , which leaves the coupling 
Hamiltonian ~H SE as a perturbation. The single-particle 
Green’s function of the disconnected Hamiltonian is de¬ 
noted by g(ui), which obviously does not mix the dis¬ 
connected regions. For the noninteracting environment, 
the respective block entries of g(ui) are available analyti- 
cally i 94 ' 116 For the interacting part the respective entries 
of g(co) are calculated via the Lehmann representation 
with respect to the reference state. This can be com¬ 
puted e.g. based on the Band Lanczos methodAii _ -ii^ 
The full steady-state Green’s function in the CPT ap¬ 
proximation is found by reintroducing the inter-cluster 
coupling perturbatively 

(5(w) -1 = <?(w) -1 — M ■ M r = M a =M,M k = 0, 

( 4 ) 

where we denote by the matrix M the single-particle 
Wannier representation of / h SE . CPT is equivalent to 
using the self-energy E of the disconnected Hamiltonian 
as an approximation to the full self-energy. Therefore, 
the quality of the approximation can in principle be sys¬ 
tematically improved by adding more and more sites of 
the leads to the central cluster. However, in doing so the 
complexity for the exact solution of the central cluster 
grows exponentially. Independent of the reference state, 
this scheme becomes exact in the noninteracting limit. 


B. Born-Markov equation for the reference state 


In the following we outline how to obtain the refer¬ 
ence state p s by using a Born-Markov-secular (BMsme), 
or more generally a Born-Markov master equation 
(BMme)^— Although this approach is standard, for 
completeness we present here the main aspects and no¬ 
tation. We loosely follow the treatment of Ref. 


The real time r evolution of the full many-body den¬ 
sity matrix p is given by the von-Neumann equation 

p = —i "H, p Typically the large size of the Hilbert 


space of R prohibits the full solution in the interact¬ 
ing case. One thus considers the weak coupling limit 
\n SE \ -C \H e \ and performs a perturbation theory in 
terms of |#?E| 80,120 

In the usual way one obtains an equation for the re¬ 
duced many-body density matrix of the system p s (r) = 
tr e {p} by working in the interaction picture pi(r) = 
e +i(y. +u ) T p(0)e -1 ^ +R ' )T with respect to the cou¬ 
pling Hamiltonian, Eq. (TXcl) . One then performs three 
standard approximations: i) Within the Born approxi¬ 
mation, valid to lowest order in \R SE \, the density ma¬ 
trix is factorized pi(r) k, /5f(r) ® pf. Furthermore, the 
environment pf is assumed to be so large that it is not 
affected by \R SE \ and thus independent of time, ii) The 
Markov approximation implies a memory-less environ¬ 
ment, that is, the system density matrix varies much 
slower in time than the decay time of the environment 
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correlation functions C a p(r). Upon transforming back to 
the Schrddinger picture this yields the BMme, which is 
time-local, preserves trace and hermiticity, and depends 
on constant coefficients, iii) To obtain an equation of 
Lindblad form which also preserves positivity one typi¬ 
cally employs the secular approximation, which averages 
over fast oscillating terms, yielding the BMsme / 8 ’ 121 ! 122 

The system-environment coupling can be quite gener¬ 
ally written in the form 'H SE = S a ®E a , with S a = S j) 

a 

and E a = E^ a . This hermitian form is convenient for fur¬ 
ther treatment.The tensor product form can be achieved 
even for fermions by a Jordan-Wigner transformation^ 
see App. [B] For our coupling Hamiltonian, Eq. (licit and 
particle number conserving systems, the coupling opera¬ 
tors take the form 


Slier — lAio- — C Aicr) (^) 

S2ia = ^=(/i<r ~ ficr)’ = ~^^ CXia ~ C Ler) ■ 

In the energy eigenbasis of the system Hamiltonian 
H s \a) = uj a |a), the BMme in the Schrodinger repre¬ 
sentation reads' 


p S (r) = -i n s + n LS ,p s (r) + £ ~ ab , cd 

1 


abed 


| a) (6| p° (r) |d> (c| - 


\d) (c| | a) (b\,p b {T) 


( 6 ) 


with 


“ ab,cd — E Ca/3 (wba,Wdc) H Sp \b) (c| S a |d)* , (7) 

a/3 


where tuba = <^b ~ w a . The Lamb-shift Hamiltonian 'H LS 
and the environment functions £ Q/ 3 (wi, LO 2 ) are defined 
in App. [Aj When employing the secular approximation, 
the terms in the BMsme simplify and in Eq. 0 one can 
replace ^apipibai^dc) ^ £,ap(bjb bj a )6 lJ j h — UJa ^ d — U j c . Hue 
to the secular approximation the BMsme can only lead to 
interference between degenerate states. The more general 
BMme also couples non-degenerate states at the cost of 
loosing the Lindblad structure of the Qme, see Sec. IIVBI 
and Ref. 0- 


where i, j denote indices of system sites. The advanced 
component follows from g A = and the Keldysh 

component g K of the finite, unperturbed system is not 
relevant for the CPT equation, Eq. 0. Once g is ob¬ 
tained, the full Green’s function is again approximately 
obtained within CPT by Eq. 0. Notice that for 17 = 0, 
G is independent of the reference state, which is why 
stsCPT, stsVCA as well as meCPT coincide (and become 
exact) in the noninteracting case. 


C. Numerical implementation 


From a numerical point of view, the two main steps are 
to first obtain the reference state p s by solving the Qme 
and then to evaluate the Green’s functions using Eq. 0 
and Eq. 0 . For the solution of the BMme, Eq. 0 one 
needs to carry out the following: i) Full diagonalization 
of the interacting system Hamiltonian which is done in 
LAPACK , making use of the block structure in N and S z . 
ii) Evaluation of the coefficients of the BMme in Eq. 0 , 
which involves coupling matrix elements {a\ S a |6) and 
numerical integration of the bath correlations functions, 
see App. |A][Cl for which an adaptive Gauss-Kronrod 
scheme is employed, iii) The steady-state p s is finally ob¬ 
tained from the unique eigenvector with eigenvalue zero 
of Eq. 0, which we determine by a sparse Arnoldi diag¬ 
onalization. Again, a block structure is related to N and 
S z . The numerical effort for the exact diagonalization 
scales with the size of the Hilbert space, and therefore 
exponentially with the system size N$- In the second 
major step, the Green’s function of the disconnected sys¬ 
tem is calculated by Eq. 0. Finally, the meCPT Green’s 
function G(oj) is found using Eq. 0 . We outline how 
to evaluate observables within meCPT and the Qme in 

App. El 


IV. RESULTS 


Single-particle Green’s function 

As discussed above, for meCPT, the Green’s function 
g(uj) of the isolated system is evaluated from the refer¬ 
ence state p s . The retarded component Eq. 0 takes the 
explicit form 

9ij(*)( UJ ) = Y.P*b x ( 8 ) 

abc 

i b \ fia |c) (cl f) G | a) + (6| f} a | c) (c| f ia I a) \ 

U! + i0 + — (u>c — <jjb) bJ + i0 + — (u>a — bJc) J 


In this section we present results obtained from the 
meCPT approach. In all calculations, except the ones in 
Sec. IIV B1 the secular approximation is applied for the 
reference state p s . The main improvements of meCPT 
with respect to bare BMsme are i) the inclusion of lead in¬ 
duced broadening effects, ii) the correct U = 0 limit and 
iii) a correction for effects missed by an improper treat¬ 
ment of quasi degenerate states in the BMsme (see be¬ 
low). In comparison to the previous “standard” stsCPT, 
meCPT also captures current blocking effects, which are 
discussed in detail in Ref. and Ref. Ho] within a Qme 
treatment. 
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left lead central molecule right lead 



FIG. 1: (Color online) Quantum dot diode: Schematic repre¬ 
sentation, see Sec. II V Al Single quantum dot with Hubbard 
interaction U and gate voltage Vg (particle-hole symmetric at 
Vg = 0), coupled via r L / R = | to a left and right lead. The 
right lead is fully polarized, i.e. only spin-f DOS is present. 
An external bias voltage Vs shifts the chemical potentials by 
Ml/r = ±^r- The leads are in the wide band limit and at 
the same temperature T. 

A. Quantum dot diode 

We first discuss a quite simple model system: a quan¬ 
tum diode based on electron-electron interaction effects. 
Fig.[I]depicts this junction consisting of a single interact¬ 
ing orbital described by a Hubbard interaction and an 
on-site term to allow for a gate voltage Vg — ^ 



where h£ = /J/ CT . The environment Eg. (llbl) . consists 
of two spin dependent, conducting leads. We model 
both, the left (L) and the right (R) lead by a flat DOS 
with local retarded single-particle Green’s function^ 

9l/r M = - 2 S ln ( ^o++j ) , with a half-bandwidth D 
much larger than all other energy scales in the model, 
mimicking a wide band limit. We keep both leads at the 
same temperature T B = T B = T and at chemical poten¬ 
tials fiL = —Hr = =r corresponding to a symmetrically 
applied bias voltage Vg. The right lead is fully spin po¬ 
larized, i.e. tunnelling of one spin species (f) into the 
right lead is prohibited while both spin species can tun¬ 
nel to the left lead. The system is coupled to the two 
leads via a single-particle hopping amplitude t' in 'H SE , 
Eq. (Hell which results in a lead broadening parameter of 

ri = ri = T^ = § = n\t | 2 ^, Eq. & and = o. 

We use F without an argument for T(u; = 0) as defined 
in Eq. (IC1I) . For meCPT we use T~L SE , see Eq. £0 , as 
perturbation. 

Such a system could be realized in: i) A “metal - artifi¬ 
cial atom - half-metallic ferromagnet nano structure 
where spin-f DOS is present at the Fermi energy while 
the respective spin-f DOS is zero, ii) A graphene nano 
structur e 31 ! 32 with ferromagnetic cobalt electrodes— iii) 
A one dimensional optical lattice of ultra cold fermions in 


a quantum simulator*^ where the hopping of spin-f par¬ 
ticles into the right reservoir is suppressed. For all three 
systems spin-f particles cannot reach the right lead, in 
the first two due to a vanishing DOS, in the third one 
due to a vanishing tunnelling amplitude. 

We consider parameters such that the junction is op¬ 
erated in a single electron transistor (SET) regime,— i.e. 
temperatures above the Kondo temperature— In this 
regime we expect an interaction induced - magnetization 
mediated blocking due to the fact that the system fills 
up with spin-f particles. On the one hand they cannot 
escape, yielding a vanishing spin-f current, and on the 
other hand they suppress the spin-f occupation, at fi¬ 
nite repulsive interaction [/, resulting also in a vanishing 
spin-f current— 

Fig. [2] (A) shows the meCPT stability diagram of the 
interacting system in the Vb — Vg plane. When applying 
a particle-hole transformation for all particles, leads and 
system, along with t' —> —t 1 we easily find the symmetry 
properties 

j(-V B ,-V g) = -j(V B , Vq) , 

(n{,)(-V B ,-V G ) = 1 - (ni)(V B , V G ) . 

From the continuity equation it is clear that only spin-f 
steady-state current can flow which limits the maximum 
current to The energies ujn of the isolated quantum 
dot can be labelled by the total particle number N and 
are for Vg = U given by ujq = 0, uq = | U and W 2 = 2 U. 
This gate voltage corresponds to the dashed line, marked 
by (X) in Fig. [2] (A). The corresponding energy differ¬ 
ences Aoi = 0.517 between the single-occupied and the 
empty dot and A 12 = 1.5 17 between double-occupied and 
single-occupied dot are associated with a further trans¬ 
port channel opening as soon as the bias V B reaches twice 
their values. The meCPT result for the current exhibits 
the well known Coulomb diamond^! close to V B = 0 and 
Vg = 0, where current is hindered because all system 
energies are far outside the transport window ±-^r, see 
Eq. (ID2I) . At Vc, = 0 a current sets in at = ±|^|, i.e. 
when transport across the system’s single particle level 
becomes allowed. The point, at which the current sets 
in, shifts with Vg linearly to higher bias voltages. This 
transition is broadened oc max(r = 0.117, T = 0.02517). 
However, not only the transport window and possible 
excitations in the system energies determine the current- 
voltage characteristics. The particular occupation of the 
system states may lead to more complicated effects, such 
as current blocking. 

Our first main result is that in contrast to stsCPT the 
blocking is correctly reproduced in meCPT. The current 
blocking is visible in Fig. [2] (A) in region (Y), see also the 
detailed data in subplot (Cl). It is asymmetric in V B 
and therefore responsible for the rectifying behaviour for 
\Vq\ > |tt|- This feature is easily understood from the 
plots of the spin resolved densities in Fig. [2] (C2). In the 
region of interest, for positive V B , (n±) = 1 which hinders 
spin-f particles from the left lead to enter the system, due 
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FIG. 2: (Color online) Quantum dot diode: (A) Stability diagram, based on the total current j = (j-y) + (j_y) as a function of 
bias voltage Vb and gate voltage Vg, obtained within meCPT. Note that (jy ) = 0. Results are depicted for T = 0.025 U and 
r = 0.1 U. (Y) marks the current blocking region. The green dashed line (X) at Vg = U indicates the parameter regime for the 
panels (B) and (C). (B) Diagonal part of the reduced density matrix pf a obtained by BMsme. (Cl) Spin-")" current j-y within 
meCPT compared to BMsme. Solid lines are for the same parameters as line (X) in panel (A). Blue dashed and solid lines for 
BMsme are indistinguishable. (C2) Spin resolved densities (n-y) and (nj.) for the same parameters as in panel (Cl), see solid 
lines in the legend. 


to the repulsive interaction U and suppresses the current. 
For negative Vg, the situation is reversed. A direct com¬ 
putation of the current in the framework of BMsme, see 
App. lD 21 also predicts the blocking, which is however not 
the case if we use stsCPT based on the zero temperature 
ground state |'f'o)s- The blocking is evident in Fig.[5](B), 
where we observe that in the blocking regime, the reduced 
density is ps = ||) (||. Independent of the value of U > 0, 
the blocking sets in at the same values of Vb in meCPT 
and BMsme. Fig. [2] (Cl) shows that within BMsme this 
regime is entered after a U independent hump in the 
current while within meCPT the hump is broader and 
weakly U dependent. The current blocking disappears 
at a bias voltage Vb oc U in both methods. Immedi¬ 
ately apparent are the much broader features in meCPT, 
which leads to a less pronounced effect in contrast to 
the total blocking predicted by BMsme. In BMsme the 
broadening parameter T enters merely as prefactor of the 
current, and broadening is solely induced by the temper¬ 
ature. This temperature induced broadening is correctly 
taken into account in both methods. For T > T the lat¬ 


ter dominates and the meCPT results are similar to the 
plain BMsme solution. A comparison of the three meth¬ 
ods is given in Tab.Q] In this simple model the blocking 
can be captured even by a straight forward steady-state 
mean-field theory in the Keldysh Green’s function with 
self-consistently determined spin densities or in stsVCA. 
This is not the case for the more elaborate system studied 
in the next section. 

TABLE I: Comparison of steady-state cluster perturbation 
theory (stsCPT), the Born-Markov-secular master equation 
(BMsme) and the quantum master equation based stsCPT 
(meCPT) with respect to their ability to capture temperature 
(T) or lead (T) induced level broadening, current blocking and 
whether the noninteracting limit is fulfilled. 


method T-broadening T-broadening blocking U = 0 


stsCPT 

yes 

yes 

no 

exact 

BMsme 

yes 

no 

yes 

approx. 

meCPT 

yes 

yes 

yes 

exact 
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B. Triple quantum dot 

In this section we discuss a more elaborate model 
system: a triple quantum dot ring junction which fea¬ 
tures negative differential conductance (NDC) based on 
electron-electron interaction effects mediated by quan¬ 
tum interference due to degenerate states as outlined in 
detail in Ref. EM . Fig. [3] (A) depicts the triple quantum 
dot ring junction, described by the following Hubbard 
Hamiltoniar422 

u s = EE”*> 

i —1 2=1 a 

+*EE/U>- ( 9 ) 

<*j> ° 

In addition to the model parameters described in 
Sec. IIV A1 a nearest-neighbour {if) hopping t is present. 
The environment, Eq. (|lbl) and coupling, Eq. (Ilcl) are now 
both symmetric in spin. Moreover, we use /j,l = —/t R = 
Y£-^T = T l = T r and T L =T R = ^ = n\t'\ 2 ^. 

Such a junction can be experimentally realized: i) Via 
local anodic oxidation (LAO) on a GaAs/AlGaAs het- 
erostructuro22 which enables tunable few electron con¬ 
trol. 30 ii) In a graphene nano structure ! 31 ’ 32 Experimen¬ 
tally the stability diagram has been explored^ alongside 
characterisation and transport measurements ! 29 ' 35 ’ 36 The 
negative differential conductance has been observed in a 
device aimed as a quantum rectifier 4- 27 Theoretically the 
study of the nonequilibrium behaviour of such a device 
has become an active field recentlyi 45 ' 46 ’ 128 4 41 

We investigate transport properties for values of the 
parameters such that the junction is in a single electron 
transistor (SET) regime, 3 ' i.e. temperatures above the 
Kondo temperature 42 In this regime we expect an inter¬ 
action induced - quantum interference mediated blocking 
as discussed in Ref. EM The rotational symmetry en¬ 
sures degenerate eigenstates labelled by a quantum num¬ 
ber of angular momentum. In situations where these de¬ 
generate states participate in the transport they provide 
two equivalent pathways through the system and lead to 
quantum interference 42 The blocking sets in at a bias 
voltage, where the degenerate states start to participate 
in the transport. It then becomes possible that a super¬ 
position is selected which forms one state with a node at 
the right lead. In the long time limit this state will be 
fully occupied while the other one will be empty due to 
Coulomb repulsion, for reasons very similar to the ones 
discussed in the previous section ! 39 ’ 40 

The steady-state charge distribution and current- 
voltage characteristics of the interacting triple quantum 
dot are presented in Fig. 0 (B, C) in a wide bias voltage 
window. The current, depicted in panel (C), in general 
increases in a stepwise manner and is fully antisymmet¬ 
ric with respect to the bias voltage direction. A block¬ 
ing effect occurs at V R ~ 1.5 |£| as can be observed in 
the BMsme and meCPT data. The previous version 
of stsCPT based on the pure zero temperature ground 


state |To )g misses this region of NDC. In contrast to 
the simpler model presented in the previous section, a 
self-consistent mean-field solution does not capture the 
blocking effects correctly in this more elaborate system. 
The BMsme solution shows many more steps in the cur¬ 
rent than the stsCPT one, which is due to transitions in 
the reference state p s of the central region. The meCPT 
results in general follow these finer steps, correcting their 
width to incorporate also lead induced broadening effects 
in addition to the pure temperature broadening. As can 
be seen in panel (Bl), meCPT predicts a large charge 
increase at the site connected to the high bias lead. Note 
that the charge density at site 2, which is connected to 
the right lead is simply: (n 2 )(Vs) = (m)(—Vg). The 
charge density at site 3 is symmetric with respect to the 
bias voltage origin. 

Next we study the impact of a gate voltage on the 
blocking. Results obtained by meCPT are depicted as 
stability diagram in Fig. [4] Upon increasing | Vg |, the on¬ 
set of the blocking shifts linearly to higher V R (Y). We 
find a Coulomb diamond for 2 Vg <; V R — |t| (D). Upon 
increasing the bias voltage out of the Coulomb diamond, 
see e.g. line (X), a current sets in but is promptly hin¬ 
dered by the blocking so that the current diminishes after 
a hump of width oc max(T, T). Interestingly this device 
could be operated as a transistor in two fundamentally 
different modes. In mode (Tl), at a source-drain voltage 
of |£| the current is on for a gate voltage of Vg = 0 
and off for Vq ~ 0.5 |t| due to the Coulomb blockade. In 
mode (T2), at a source-drain voltage of « 1.5 |£| the cur¬ 
rent is off for a gate voltage of Vg = 0 due to quantum 
interference mediated blocking and on for Vq = 0.25 \t\. 

Next we discuss the current characteristics in the vicin¬ 
ity of the blocking in more detail, as well as the impact 
of the interaction strength U. The first row of Fig. [5] 
shows the total current through the device for different 
values of U. The blocking region shifts to lower bias volt¬ 
ages with increasing U. As discussed earlier, structures 
in the BMsme results are only broadened by tempera¬ 
ture effects in the steady-state density (compare e.g. the 
width of the structures in the local density in the second 
row of Fig. ©, while meCPT additionally takes into ac¬ 
count the finite life time of the quasi particles due to the 
coupling to the leads, given by 1/r. This can be seen 
by solving Eq. © for the local Green’s function at de¬ 
vice sites. Especially for higher lead broadening T this 
gives rise to significant differences in the meCPT results 
compared to the BMsme data. From the bottom row 
of Fig. [5] we see that, before the blocking regime is en¬ 
tered, the steady-state changes from a pure N = 2 state 
to a mixed N = 2 / N = 3 state at the hump in the 
current. Obviously, blocking arises because the system 
reaches a pure N = 3 state for U = 2 \t\ and U = 3 |t| 
at V R s=s 1.4 |t|. For U = \t\ the current is only partially 
blocked, because the contribution of the N = 2 state is 
not fully suppressed. For all [/-values, however, we find 
NDC. As far as the meCPT current is concerned, the 
complete blocking at higher interaction strengths, pre- 










FIG. 3: (Color online) Triple quantum dot (A) Schematic representation, see Sec. IIV Bl System Hamiltonian as defined in 
Eq. ([9]). Site 1 couples to the left lead and site 2 to the right one, both with r L / R = The leads are held at the same 
temperature 7 l/r = T and the chemical potentials /x L / R = are shifted by the bias voltage. (B) Local charge density (ni) 
as a function of bias voltage Vb ■ The results are obtained by meCPT, BMsme and stsCPT, see color code of panel (C). (C) 
Total current j = ^2 a {jLio-) into the system at site 1 as a function of bias voltage Vb- Results, shown in panels (B,C), are for 
U = 2 |t|, T « 0.02 |t|, r = 0.1 |f| and Vg = 0, corresponding to line (X) in Fig. [4] 



- 3 - 2-10123 



FIG. 4: (Color online) Triple quantum dot: stability diagram. 
Total current entering the system as a function of bias voltage 
Vb and gate voltage Vg, obtained within meCPT. The block¬ 
ing region is indicated by (Y), the Coulomb diamond by (D). 
The two arrows (Tl) and (T2) mark two device operation 
modes as discussed in the text. All results are for U = 2 |£|, 
T = 0.02 |f| and F = 0.1 |f|. Dashed line (X) for Vg = 0 marks 
the parameter region depicted in Fig.QT] (C). 


dieted by BMsme, is reduced to a partial blocking due to 
the lead induced broadening effects in meCPT. Although 
p^ b changes significantly twice in the blocking region (for 
U = 2 and U = 3), the charge density (rij) just increases 
once from (m) ~ 0.75 to (ni) « 1. 

Details of the steady-state dynamics are provided in 
Fig.m Before the blocking region is entered (Vb = 0.4 |i|) 
the system is in a pure state with N = 2, which corre¬ 
sponds to the zero temperature ground state | v E'o),s i n 
the N = 2 sector. Here the transmission function T(w), 
Eq. (ID3I) . of meCPT agrees with the one of stsCPT. A 


small current is obtained due to the N = 2 —> 3 excita¬ 
tion at uj ss 0.55 |f|. Increasing the bias voltage has no in¬ 
fluence on the reference state in stsCPT, which therefore 
remains in the N = 2 particle sector. Consequently, the 
transmission function in stsCPT does not change. Only 
the transport window increases linearly with increasing 
Vb- For Vb = 1.4 \t\ it includes the peak at ~ 0.7|f| and 
results in a significant increase in the current obtained 
in stsCPT (see stsCPT result in Fig. [3]). This is in stark 
contrast to the BMsme current, depicted in Fig.[5j which 
exhibits perfect blocking for Vb = 1.4 |i|. The reason 
for the current-blocking is that only two states, both in 
the N = 3 sector and doubly degenerate, have significant 
weight in p^ b . The meCPT solution is based on the mod¬ 
ified density matrix and therefore the current is dimin¬ 
ished, since the next possible excitation is at w ~ 0.9 |t| 
(.N = 2—^3), which is outside the transport window 
W(uj) « (—0.7|£|, 0.7|f|), Eq. (ID2I) . Due to the lead in¬ 
duced broadening of T(w) and the temperature induced 
broadening of the transport window, the current is how¬ 
ever only partially blocked. For Vb = 2.4 \t\ this excita¬ 
tion falls into the transport window and the current is 
no longer blocked. In this case, the state is a mix¬ 
ture of N = 2, 3,4. The dominant excitation responsible 
for this current is again the ground state excitation at 
u> ~ 0.55 |t| from N = 2 —> 3. This is why in this regime 
the stsCPT current, based on the pure two particle state 
is again similar to the meCPT current. 

Our results on the Qme level have been checked with 
those presented by Begemann et al. in Ref. [3£j and Darau 
et al. in Ref. [13 for a six orbital ring which shows similar 
blocking effects. Different types of blocking effects in 
various parameter regimes have been discussed in detail 
in a Qme framework also for the three orbital ring by 
Donarini et al. in Ref. M 
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U = l|t| 



0 0.5 1 1.5 2 2.5 


U = 2|t| 



0 0.5 1 1.5 2 2.5 


U = 3111 



0 0.5 1 1.5 2 2.5 



FIG. 5: (Color online) Triple quantum dot: Dependence of the current blocking on the interaction strength U. (Top row) 
Total current j as a function of bias voltage Vb- (Middle row) Charge density (m) at site 1. The color code of the top row 
is valid. (Bottom row) Summed diagonal elements of the density matrix wn = X^eiv Paa P er particle number N. The black 
markers in the mid panel (U = 2 |f|) indicate for which Vb detailed results are given in Fig. [6] Solid lines in all panels are for 
T = 0.02 |f|, T = 0.1 |t| and Va = 0. Results for T = 0.1 |f| are depicted in the central panels by dotted lines and those for 
r = 0.5 |f| in the right panels by dashed lines. 


V B =0.4|t| 


1 

0.8 
0.6 
0.4 
0.2 
0 

-2.5 0 2.5 

«V|t| 

FIG. 6: (Color online) Triple quantum dot: Dynamic transmission function T(cu), Eq. (iDl . as obtained by meCPT and stsCPT. 
Same parameters as in Fig. [S] (bottom mid) at the three indicated bias voltages: Vb = 0.4 |£| (left), Vb = 1.4 |t| (middle) and 
Vb = 2.4 |£| (right). The temperature broadened transport window W(u>), Eq. (1D2I) . is depicted as a dashed black line. 
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Quasi-degenerate states 

Next we study the reliability of the secular approxima¬ 
tion in the case of quasi degeneracy of the isolated en¬ 
ergies of the system and benchmark its applicability to 


create a reference state for meCPT. To this end we ap¬ 
ply a second gate voltage that couples only to the third 
orbital, see Fig. [3] (left), and leads to an additional term 
Vg ,3 in the system Hamiltonian. This lifts the degen¬ 
eracy of states present at Vg ,3 = 0 and therefore requires 
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a treatment within the BMme, see Ref. Ho|. 

In the following we discuss the same parameter regime 
as above. In Fig.Q we present results obtained using 
meCPT (solid lines) and Qme results (dashed lines) for 
the BMsme (A) and for the BMme (B). The meCPT re¬ 
sults of each panel are obtained using the respective Qme. 
In the BMsme data a very small |Vo, 3 1 has a drastic ef¬ 
fect on the current-voltage characteristics. The block¬ 
ing present at Vq ,3 = 0 is immediately lifted by very 
small |Vb,s| and the current jumps to a plateau. For 
larger |Vg, 3 | the current stays on this plateau until fur¬ 
ther transport channels open up. This ”jump“ at small 
|Vg, 3 | arises due to the improper treatment of quasi de¬ 
generacies in BMsme. MeCPT results based on BMsme 
show a smooth change of the current-voltage character¬ 
istics. BMme on the other hand correctly accounts for 
the coupling of the quasi-degenerate states and also ex¬ 
hibits a smooth dependence on Vg, 3 . For meCPT based 
on BMme we find qualitative similar results to meCPT 
based on BMsme, which emphasizes the robustness of 
the meCPT results in general. From this it is appar¬ 
ent that meCPT is capable of repairing the decoupling 
of quasi-degenerate states in the BMsme to some degree. 
However, to study blocking effects at quasi degenerate 
points it is of advantage to make use of the BMme in 
meCPT. 

As discussed below in Sec. II V Cl the BMme is not of 
Lindblad form and does not necessarily result in a pos¬ 
itive definite reduced many-body density matrix p^ b in 
general. Using a not proper density matrix in Eq. © 
may result in non-causal Green’s functions when the 
steady-state p^ b is obtained from the BMme. This can 
be avoided by using a modified reference state —> 

pf fc 0(A — |w a — W(,|), with 0(x) the Heaviside step func¬ 
tion and A a small quantity, being e.g. ss 10 -6 , in Eq. ([5]). 
which renders the Green’s functions causal. This is some¬ 
what an ad-hoc approximation and should be seen simply 
as a way to explore the effects of continuously breaking 
degeneracy in the problem. 


C. Current conservation 

Finally we comment on conservation laws in meCPT. 
Within BMsme and BMme the current conservation 
(continuity equation) is always maximally violated in a 
sense that the current within the system is zero. This 
is due to the zeroth order p s as discussed in App. ID 21 
In BMsme the inflow from the left lead into the system 
however always equals the outflow from the system to the 
right lead. Without the secular approximation the quan¬ 
tum master equation (BMme) is not of Lindblad form 
and the final many-body density matrix is not guaran¬ 
teed to be positive definite. 121 ! 122 This in turn can lead 
to slightly negative currents in regions where they are 
required to be positive by the direction of the bias volt- 
ageS. Furthermore, the inflow can be slightly different 
from the outflow. 


(A) Born-Markov-secular 



v B /|t| 


(B) Born-Markov 



v B /|‘l 


FIG. 7: (Color online) Triple quantum dot: Effects of lifting 
degeneracies in the system energies by a third gate voltage. 
Total current j as a function of bias voltage Vb , for three 
different gate voltages Vg ,3 applied to site 3. Results based on 
the Born-Markov-secular approximation are compared with 
those of the Born-Markov approximation. Solid/dashed lines 
indicate the meCPT/BM(s)me result. All results are for U = 
3 |£|, T = 0.02 |t| and F = 0.1 \t\. 


In the noninteracting case, meCPT fully repairs the vi¬ 
olation of the continuity equation present in the reference 
state. For increasing interaction strength, the violation 
of the continuity equation typically grows also in meCPT. 
In particular, the overall symmetry of the current stays 
intact (in our case, inflow equals outflow), while the cur¬ 
rent on bonds between interacting sites does not exactly 
match the current between noninteracting sites. This 
typically small violation of the continuity equation can 
be attributed to the violation of Ward identitie s 135 ’ 136 in 
the non-conserving approximation scheme of CPT. 137 ' 138 


V. SUMMARY AND CONCLUSIONS 

We improved steady-state cluster perturbation theory 
with an appropriate, consistent reference state. This ref¬ 
erence state is obtained by the reduced many-body den¬ 
sity matrix in the steady-state obtained from a quantum 
master equation. The resulting hybrid method inher¬ 
its beneficial aspects of steady-state cluster perturbation 
theory as well as from the quantum master equation. 


























11 


We benchmarked the new method on two experimen¬ 
tally realizable systems: a quantum diode and a triple 
quantum dot ring, which both feature negative differen¬ 
tial conductance and interaction induced current block¬ 
ing effects. meCPT is able to improve the bare quantum 
master equation results by a correct inclusion of lead in¬ 
duced level-broadening effects, and the correct noninter¬ 
acting limit. In contrast to previous realizations of the 
steady-state cluster perturbation theory, meCPT is able 
to correctly predict interaction induced current blocking 
effects. It is well known that the secular approximation 
(BMsme) is not applicable to quasi degenerate problems, 
which is corroborated by our results for the steady-state 
current. However, meCPT based on the BMsme density, 
is able to repair most of the shortcomings of BMsme. The 
results are very close to those obtained by meCPT based 
on the density of BMme, where the quasi-degenerate 
states are treated consistently. 

The computational effort of meCPT beyond that of 
the bare quantum master equation scales with the num¬ 
ber of significant entries in the reference state density 
matrix but is typically small. In the presented formu¬ 
lation the new method is flexible and fast and therefore 
well suited to study nano structures, molecular junctions 
or heterostructures also starting from an ab-inito calcu- 
lationJ^ 
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poo pO 
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J 0 J — oo 

(A3) 
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OO 
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—OO 

oo oo 

A«/ 3 (w) = [ dr sign(r )C a p (r)e +lulT = - if dw' ^ af) ^] . 

J 7 T J U) -U)' 

—oo —OO 
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Appendix A: Born-Markov and Pauli master 
equation 


Here we provide the detailed expressions for the co¬ 
efficients in the BMme and BMsme of Eq. 13 and dis¬ 
cuss the equations governing the time evolution into the 
steady-state. 

The Lamb-shift Hamiltonian is defined as H LS = 
J2 A ab | a) {b\, with 

ab 

A-ab — ^ ^ ^ ^oipi^bci ^ac) (c| |^) (c| S& |tt) 

0-/3 c 

(Al) 


Note that [H LS ,i-L s ]- = 0. In the secular approxima¬ 
tion (BMsme) one can replace A Q t p(u)bc,u ac ) —> A a p{pJ b — 
0J c )bu> b ,uj a ■ The expressions for the BMme and BMsme 


Eq. 13 are valid if 


n E ,p E 


= 0 and tr j = 0. 


The coupled equations for the real time evolution of the 
components of the reduced system many-body density 
matrix p$ b = (a| p s \b) according to the BMsme read 

Pab( T ) = - Wa)pfb(r) (A7) 

+ (pac( T ) A cb - KcP S cb (r)j 

c ^ ' 

d~ 'y 1 ( i= ‘ac,bdPcd{ T ) ~ ^cd,caPdb ( T ) 

cd \ 

~ 2 ^cb,cdP a d{ T ) 

The equations simplify further for system Hamiltoni¬ 
ans 7 I s with non-degenerate eigenenergies u) a . Then the 
diagonal components (j> a = pf a decouple from the off- 
diagonals and one recovers the Pauli master equation for 
classical probabilities 

< M r ) = f S “^ c (r) - Eclair) J , (A8) 

C ^ ' 

with simplified coefficients 
^ab * = “ ab,ab = ^ ^ ^a) (&| Sft |6) (tt| S a \b ) 

a>/3 
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In this case the dynamics of the decoupled off-diagonal 
terms (a ^ b) is given by 

pfb(r) = + A b-Wa — A a ) 



where the simplified Lamb shift terms are 

A a • — A aa — ^ ^ ^ a ^c) (c| Sfi |u) (c| S a |ct) 

ck/3 c 


Appendix B: Hermitian tensor product form of the 
coupling Hamiltonian 

For the BMsme (see Sec. IIII Bl) it is necessary to bring 
the fermionic system-environment coupling Hamiltonian, 
Eq. (licit to a hermitian tensor product form, which re¬ 
quires [S a ,E a ]_ = 0. For the fermionic operators in 
Eq. (O we however have [/Lc Afc(T ]^ = 2 fl a c xkr7 . A so¬ 
lution is provided in Ref. l79j by performing a Jordan- 
Wigner transformation^^ on the system and environ¬ 
ment operators 

fia = ]^[ (£l ® ■ • ■ <8> £|-l£i M+1 ® ■ ■ ■ ® IIls)^^ 
a 

® (^1 ® ■ ■ ■ ® ^ L e )e,\<t ) 

A 

<7 

II (rj l ® ... ® 1 ® • ■ ■ ® , 

A 

where £* and T)j denote local spin-i degrees of freedom 
at the system and environment sites respectively and the 
overall ordering of operators is important. Ls/Le denote 
the size of the system / environment. Reordering Eq. © 
we find nf E = J2 t' Xij< 7 flc Xj<7 - where the 

ijcr 

minus sign arises due to the fermionic anti-commutator. 
Plugging in the Jordan-Wigner transformed operators 
leads to 


where in the last line we have defined new operators 

fia = [£j ® [~fi +i ® • • • ® £l s ]] a > 

c\ia = *A ija IK ® • • ■ ® Vj-l] ® vy ] Xa . 

3 

A* = 5Z [vf ® [»?* ® • • ■ ® ^-i]] Act ■ 

3 

Note that the phase operator Pi(jX)<r = 

[-£?+i ® • • • ® ££ s ® faf ® • • • ® ^I-i]a] 

!+E E n m +N jX i 

(—1) x ' m=i + 1 counts the particles between 

system site i and environment site j for spin a de¬ 
pending on the ordering of the environments A. It 
is straight forward to show that the bar operators 
fulfil fermionic anti-commutation rules. Furthermore 
[fia,c X i a ]- = 0, which allows us to write the coupling 
Hamiltonian in a tensor product form. Note that in 
general [fia,cyja]- ^ 0 for i ^ j which is however not 
relevant for the coupling Hamiltonian where only the 
same i couple. 

The new operators in hermitian form are given in 
Eq. (J5]) by replacing c —> c and /—>■/. Next we show, by 
examining the BMsme, that in most cases the additional 
phase operator in c drops out of the calculations and we 
are even allowed to use the original / and c operators 
instead of the barred ones. The operators c only enter 
the equations in the environment correlation functions 
C a p{j) as defined in Eq. (|A2j) . Plugging in the barred 
operators we obtain for normal systems which preserve 
particle number 

C aP {T) oc tr {e +mET fl j<7 e- mET P? {jX)<T c Xj<7 p E } , 

with = 1, where we required that [TL E ,Pij]- = 0. 
The dropping out of the phase operators implies that 
for normal systems where the disconnected environments 
conserve particle number we can omit the Jordan-Wigner 
transformation and do all calculations as is with the orig¬ 
inal environment creation/annihilation operators in her¬ 
mitian form. 


Appendix C: Bath correlation functions 




IJ(J 


'Xija 




[-£f+i ® .. • ® ££ s ® Vi ® ■ ■ • ® Vj-i] ® Vf 

c ® [-er+i ® • ■ • ® ?£ s ® ® • 


+ t 


- a X 


Xija 


Vj-l) ® vt 


E (a 


- a X 


^ ** Ci + fi ® 4 


In the wide band limit, analytical expressions for the 
bath correlation functions are available in Ref. [39|. For 
arbitrary environment DOS, explicit evaluation of the en¬ 
vironment correlation functions becomes convenient for 
hermitian couplings, Eq. ([5]) as outlined in Ann. iBl 'A Es¬ 
sentially the environment functions can all be obtained 
via integrals of the environment DOS p(u). Care has to 
be taken when going to very low temperatures and solv¬ 
ing the integrals with finite precision arithmetic to avoid 
underflow errors. 
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The time dependent environment correlation functions 
C al 3 ( 7 -), Eq. (1A2I) become 

OO 

Cu(t) = G 22 (t) = f^Y f dvT *a{v) 

*"-oo 

x ^e~ lVT + 2ip FF) (v, T\, p x ) sin (ut) 

OO 

Ci 2 (r) = -G 21 (t) = ^Y f duT ^i v ) 

* ^-00 

x ^ - e.~ WT + 2pfd0, T x , h x ) cos (vr)j , 
where C a p(r) = G^ a (—r) and the coefhcient 

T Xa (v) =2Tr\t' Xa \ 2 Y S ( u - “\ka), (Cl) 

k 

is proportional to the lead DOS. 

For the BMsme, the respective full even Fourier trans¬ 
forms t, a p(oj), Eq. (IA5I) we find 

£11M = £22 M = 

^Y Tx °(- UJ ')P FD (~ u ’ V\) + F Xc r(u)p FD (uj,T x , g x ) , 

A ct 


tarded CPT Green’s function 

OO / 

da i f dui 


£12 (w) — — £21 (w) — 

1 

2 


Y t >«?(-u)pfd(-u, fi x , [i x ) - T Xa (uj)p FD {uj,T X: g x ), 


Act 


where p FD (ui,T, p) = 1 -p FF) (w,T, p). 

The odd Fourier transforms A a p(co), Eq. (IA6I) are given 
by 


K ija — 2 2 / 27T l Pincr(p>)(Gj na {u)) 


+ <?£» ([-Pcr(w), M a \_) nm (Gf ma (uj))*J , 

where M a is the inter-cluster perturbation defined in 
Eq. (U). Here we use the Einstein summation conven¬ 
tion, the last line holds within CPT and Pij a { w) = 

di j (1 2p F D (iO j Ti , . 

From the real part of the single-particle density-matrix 
we read off the site occupation (rii) = K Ha the spin 

CT 

resolved occupations (ni„) = Ku a and the magnetization 
(mi) = ^(Kuf - Kui). 

The current (juj)) between nearest-neighbour sites (ij) 
is related to the imaginary part of Kjj a and reads in sym¬ 
metrized form 

g 

(j(ij) ) ~0f\ (^'b cr ^'b (T d jiaKjia) 5 

which is of Meir-Wingreen form22, and hij a is the single¬ 
particle Hamiltonian. 

Equivalently, the transmission current between two en¬ 
vironments A = 1, 2 can be evaluated in the Landauer- 
Biittiker f orm 55 1 65 1 98 

OO 

0'i/2) = | J ^W(cv)tr{T(u)} , (Dl) 

—OO 

with the transport window 

W(w) = Pfd(w, Ti, fii) - Pfd{u,T 2 , H 2 ), (D2) 


An(w) — A 22 (w) — 

OO 

i-Y d " r ^( i 

\ _ ■J 


Y Pfd(g P\,p\) _ Pfd(c P\,P\) 


v + ui 


V — OJ 


Ai 2 (w) — — A 2 i(w) — 


1 

27T 


Y if dvT *°(y) 

Act „ 


PfdYYajMa) Pfd( z/ Ya>Ma) 


v + u> 


and where the transmission function 

7» = G fl (w)ri(w) (G iJ ( w )) t r 2 (w), (D3) 

is given in terms of G R (u>) = 

^(^^(w)) -1 — (Ei + S 2 )j with the lead broaden¬ 
ing functions of lead A projected onto the system sites 
i,j is EAij = M iX g^ x M Xj and F x = -23m (e a ), 
compare also Eq. CH). 


Appendix D: Evaluation of steady-state observables 

1. Steady-state cluster perturbation theory 

Within meCPT single-particle observables are avail¬ 
able by integration of G(w), Eq. d4j). Its easy to show 
that the single-particle density matrix Kij a = -j- — 

OO 

1 / can be expressed in terms of the re- 

— OO 


2. Born-Markov master equation 

Within the Qme, basic single-particle observables are 
available in terms of the reduced system many-body den¬ 
sity matrix p s . The single-particle density matrix n reads 

Kija = tl[flf ja pA = Y ( b \ fLfja l«) Pab , ( D4 ) 

' ' ab 
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where a and b denote eigenstates of the system Hamilto¬ 
nian. Note that within the BMme/BMsme K l]a is purely 
real and therefore does predict zero current. 

However, an expression for the current to reservoir A 
can be found by making use of the operator of total sys¬ 
tem charge Q and total system particle number N, where 
q denotes the charge of one carrier 


and for non-degenerate systems, in the Pauli limit we find 
from the BMsme 


EjV) = ^QW)=« tr (^ S (r)) • 


Jnon-deg = 9 ( n “ “ ^ 

ab 


Taking p s (r) from the Qme we obtain 


j X = <lJ2 l nc ~ ) Pab > 

abc 
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